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Second-harmonic generation (SHG) from periodic arrays of subwavelength rectangular air hole with 
various aspect ratios perforated in gold thin films can get resonantly enhanced for some specific geometric 
shapes. Here we clarify the physical origin of this shape resonance effect. A nonlinear coupled-mode theory 
is set up to solve energy conversion from fundamental wave (FW) mode to second-harmonic wave (SHW) 
mode within the nanoscale air hole. It reveals that several physical mechanisms, including the FW mode 
excitation amplitude, FW-SHW modal spatial overlap, FW-SHW mode phase mismatch, and SHW mode 
attenuation, are all geometric shape sensitive and altogether act to induce the SHG shape resonance effect. 
The theory agrees well with experimental observations and provides an accurate and complete explanation 
for the long-emphasized but elusive shape effect. The study may stimulate deeper insights to visualize 
general nonlinear nanophotonic processes and pave the way to engineering high-efficiency nonlinear 
nanophotonic devices. 



Light propagation through subwavelength structures has exhibited a variety of interesting properties that are 
otherwise not observable in macroscopic materials and structures 113 . Experiments and theories have shown 
that the linear transmission intensities through sub-wavelength air-hole arrays perforated in metal thin films 
are sensitive to many factors, such as array periodicity, hole shape 811 , hole size 14 , incidence angle 15 , and film 
thickness 1617 . The physical mechanisms behind these linear- optical shape effect have been well established, where 
the excitation, transport, and coupling of surface plasmon polaritons have played the key role. 

In addition to the fascinating properties of linear light transport, nonlinear optical properties, in particular, the 
second-harmonic generation (SHG) from metallic nanostructures has also been investigated extensively, such as 
those in single nanoparticle 18 , bowtie nano-apertures and nano-antennaes 19,20 , L-shape arrays 21 , and the silver film 
with triangular air holes arrays 22 . Methods of enhancing the SHG were also studied widely, such as those by 
double plasmonic resonance 23 , Fano resonance 24 , and even nanocup structure 25 . Asymmetry and circular dichro- 
ism in optical SHG emission of gold nanostructures have also been studied 26 2S . Recently, it was further shown 
that geometric shape factors such as spacing between the double holes 29,30 and the aspect ratio (AR) of rectangular 
holes 31 critically determine the resonant enhancement of second harmonic generation (SHG) in these metallic 
nanostructures 29 33 . For instance, Ref. 31 reported an experimental observation that the SHG signal could vary by 
several orders of magnitude when the AR of rectangular holes was tuned and a remarkable resonance emerged at a 
specific value of AR. These studies clearly showed the essential role of shape effect on nonlinear optical interaction 
in metallic nanostructures. However, the fundamental mechanism of the resonantly enhanced SHG signal is so far 
not clear, as the physics is more complicated than in the linear-optical regime. To understand why shape effect 
occurs, slow light effect has been suggested in Ref. 31. However, further investigation showed that the SHG signal 
does not enhance in proportion to the reduction of the group velocity 34 . Therefore a clarified microscopic picture 
about the shape resonance effect is still lack and highly desired. 

In this paper we investigate the explicit origin of shape effect in these nonlinear optical nanostructures both 
experimentally and theoretically. In our experiment, we measured the SHG from a periodic array of subwave- 
length rectangular air holes using femtosecond laser pulses and observed the dependence of SHG signal intensity 
on the AR of the holes. Particularly we measured the propagation time duration of the pulses through the air hole 
with an aim to clarify the slow light effect. In our theory, we present a nonlinear coupled-mode theory analysis 
based on Maxwell's equations that describes the interaction between fundamental wave (FW) and the second 
harmonic wave (SHW) within the nanoscale air holes. The theory shows that the transverse modal coupling 
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and conversion between FW and SHW modes and their transport 
characteristics play an essential role in the shape-resonant effect of 
SHG observed in our experiment and previous reports 31 . Our theory 
is of great help to draw a clear picture about how several physical 
effects, including the modal profile similarity, the local field enhance- 
ment, the phase mismatch, and modal attenuation, cooperate and 
compete in ultimately shaping the SHG features in these nonlinear 
metallic nanostructures. 
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Figure 1 | Experiment on SHG from nanostructured gold thin film. 

(a) The schematic experimental setup. FM: flip mirror; CP: chopper; BBO: 
nonlinear crystal; fl and f2: focusing lens, (b) The scanning electron 
microscope image of a typical sample of the rectangular hole array, 
whereby AR = 2. The AR of the rectangular hole is defined as the ratio 
between the y-axis edge length (260 nm here) and the x-axis edge length 
(130 nm here), (c) A high resolution image, with the x-y direction and the 
laser polarization (parallel to the y-axis) marked on the figure, (d) The 
schematic physical picture illustrating generation and transport of SHW 
mode E 2 fi(x,y) from surface gold atoms of the thin film under excitation of 
the FW mode Eij)(x,y) within the air hole (which can be seen as a channel). 



Results 

The experimental setup for SHG measurement is depicted in 
Fig. 1(a). In our experiments, a series of rectangular hole array sam- 
ples of different AR were prepared in a gold thin film of 160 nm 
thickness, where the 20 X 20 square-lattice array composes of rect- 
angular holes with a spatial periodicity of 410 nm and a total area of 
8 um X 8 um (Fig. 1) (see Method, Sample Preparation). All the 
samples are designed to have the same area as S = 260 nm X 
130 nm. The scanning electron microscopy (SEM) picture of a typ- 
ical sample (with AR = 2) is displayed in Figs. 1 (b) and 1 (c). The light 
spot on the sample had a diameter of about 30 urn. As depicted in 
Fig. 1(d), under the excitation of high-peak power FW laser pulse, the 
surface of gold thin film would generate SHW, which will propagate 
through the air hole and eventually radiate out from the thin film. 
When needed, each sample was rotated by 90° around the z axis to 
facilitate taking two AR-dependent data. The measurements of SHG 
waves (see Method, Measurement of SHG signals) and of time delay 
can be easily switched. 

In Fig. 2(a) we show the observed SHG signal intensity as a func- 
tion of the aspect ratio. A prominent peak is observed for an AR value 
around 2. In the inset we show a typical spectrum of the SHG signal 
for the AR = 2.037 sample. Our data are consistent with those 
reported in Ref. 31. The time delay in group velocity was measured 
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Figure 2 | Experimental data of SHG from different gold thin film 
samples, (a) The measured SHG signal intensity versus the aspect ratio of 
the rectangular hole. The inset is a typical spectrum of the SHG signal. A 
resonance peak is clearly observed for the AR = 2.037 sample, 
(b) Convolution temporal profile of the light pulse after propagating 
through the AR = 2.037 and AR = 0.491 patterns, with a relative time- 
delay readout of 6.5 fs. The inset is the measured time delay versus aspect 
ratio, obtained by taking square vacant areas as references. 
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by cross-correlating the transmitted fundamental pulse with a ref- 
erence pulse in a 1 mm-thick beta barium borate (BBO) crystal 35 . We 
averaged for 15 scans and the overall precision obtained became 
exceptionally high, which was about 2 fs (see Methods, Method to 
achieve ultrahigh temporal resolution). 

To determine the absolute time delay, we compared the temporal 
traces through patterns with that through a large square hole on the 
Au film. A 8 X 8 um 2 vacant area were milled on each Au film. In 
such a way additional effect caused by the spatial pulse tail is mini- 
mized. The results are shown in the inset of Fig. 2(b), where we have 
subtracted the peak positions for the patterns by those for the vacant 
areas, respectively. Note that the two data points around AR = 2 are 
for two different samples. For all the AR samples discernable time 
delay was clearly observed for the laser pulses after traversing 
through the metal nanohole array, as compared with those through 
the vacant window, indicating the existence of light slowdown. 
However, overall there does not exist a dramatic slowdown for any 
of the samples. Comparison of the value for AR = 0.5 with that for 
AR = 2.0 best tests the slow light effect, since most of the other 
factors such as the laser tail effect are absent for a rotated sample. 
Also for this reason, to elucidate a pure shape effect, it is better to 
compare the AR = 2 and AR = 0.5 cases, rather than to compare the 
AR = 2 case with air. 

In Fig. 2(b) we show explicitly the temporal traces for the AR = 
2.037 and AR = 0.491 patterns, respectively. Note that there is only a 
6.5 fs delay between these two cases. The delay for the AR = 2.037 
sample is 2 ~ 3 times larger than that for the AR = 0.491 pattern [see 
inset of Fig. 2(b)]. However the bit rates of slowing-down in group 
velocity, Afd e i ay /Afp U i se( j uration , are 0.065 and 0.208 respectively, both 
of which are much less than 1. Therefore the nonlinear interaction 
time for the AR = 2.037 sample is not dramatically longer than that 
for the AR = 0.491 sample. Apparently the slowing-down exists 
[inset of Fig. 2(b)], however, it does not appear to be the major factor 
leading to the shape-dependent resonance. An alternative explana- 
tion for the fundamental mechanism of the nonlinear shape res- 
onance effect must be found. 

Maxwell's equations, being scale invariant, provide a secure basis 
for optical problems ranging from classical optics to fashioned nano- 
photonics, from linear optics to nonlinear optics. Therefore, they are 
naturally the starting point for accurately describing various non- 
linear nanophotonic problems and surely will provide the ultimate 
answer to our puzzling problems of shape resonant SHG from metal- 
lic nanostructures. In conventional nonlinear optics for macroscopic 
materials two factors are crucial for generating nonlinear optical 
signals: energy conservation and phase matching (thus momentum 
conservation). However, for nanostructures, there might be other 
crucial factors that dictate the nonlinear optical processes but have 
been overlooked in conventional nonlinear optics. 

Before we go into the details of the theoretical analysis, let us first 
draw a general physical picture about the SHG from the metallic 
nanostructure, which has been depicted in Fig. 1(d). The first thing 
is the source of SHG. SHG comes from a second-order nonlinear 
process originated from the surface atoms of nanostructured gold 
thin film, where the broken spatial inversion symmetry will lead to 
non-zero second-order nonlinear susceptibility. The illumination of 
high peak power femtosecond pulse laser upon these surface gold 
atoms will induce second-order nonlinear polarization, which will 
radiate SHW. The second thing is that in principle all the surfaces, 
including the outer surfaces (the front and back side of gold film) and 
the inner surfaces (the wall of each rectangular air hole), will con- 
tribute to SHG. However, as the linear electric field will be subject to a 
much larger enhancement effect within cavity-like nanoscale air 
holes than in the front and back side of the gold film, it is expected 
that the SHG within the air holes will make a dominant contribution 
to observed SHG signals, and this is the first stage (the generation 
process of SHW) where the geometric shape of air hole can come in 



to play a role. The third thing is that the radiated SHG signals must be 
subject to the modulation by the air hole before they can ultimately 
transport and transmit through the nanostructure gold thin film. 
This is the second stage (the transport process of SHW) where the 
geometric shape of air hole comes in to play a role. At this stage, both 
the transverse modal profile and longitudinal modal dispersion of the 
F W and SHW within the air hole can make significant contributions. 
The final thing is that the detected SHG signal is just the summed 
contribution from all the SHW exiting the air hole and enduring 
zero-order Bragg diffraction. The above analysis has decomposed 
the whole complicated nanoscale SHG problem into several pro- 
cesses with the corresponding underlying physical mechanism being 
easy to understand. By this, the problem becomes much simpler to 
yield a much clarified picture, allowing the major physical mech- 
anism of the puzzled shape resonance effect to be disclosed. 
However, the ultimate answer to this question again relies on quant- 
itative solution and analysis based on rigorous Maxwell's equations. 

According to the above argument, the SHG within each air hole 
may play the key role. Besides, as adjacent air holes do not have direct 
optical coupling due to strong attenuation of light within bulk gold 
materials, we can consider SHG within a single air hole. By the way, it 
is worthy to mention that the film thickness (160 nm) is comparable 
with that of the lateral dimension (260 nm X 130 nm) of the metallic 
hole. Thus it is safe to take the hole as a channel or waveguide for the 
light propagation. For these reasons, we start by applying Maxwell's 
equations to a single unit air hole taken as a waveguide (although 
with a very short length). It is assumed that the rectangular air hole 
has its edges parallel to the x- and y-axis, respectively, and the light 
transports along the z-axis. The transport of linear light through the 
air holes and the corresponding three-dimensional (3D) electric field 
distribution can be readily solved by 3D finite-difference time- 
domain (FDTD) method (see Method, Numerical simulations). 
The wave equation for the SHW field within the air hole can be 
written as 
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term of SHW, and Ei and E 2 denote FW and SHW fields, respect- 
ively. The nonlinear susceptibility y} 2 ' 1 is a second-rank tensor in 
general, which is non-zero only at the surfaces where spatial sym- 
metry is broken. As the outer and inner surfaces of our gold nanos- 
tructure samples have an amorphous morphology (see Method, 
Sample preparation), the general second-rank tensor y} 2 ' 1 reduces 
to an isotropic quantity that contains only non-zero diagonal ele- 
ments T 2 } r = yj-2, = y£L = y„, which is assumed to be uniform for all 
'^xxx yyy '^zzz 

the surfaces. Thus we 

haveP (2) = x!, 2) |Ei! 2 - For 

none-depletion non- 
linear optical processes, this can be seen as a perturbation term in the 
inhomogeneous differential equation Eq. (1). 

In our solution we only consider the contribution from the dom- 
inant metallic-waveguide mode (with the smallest attenuation 
strength) for both FW and SHW fields. Other higher order modes 
have much larger attenuation strengths and are expected to make 
negligible influences to the solutions. With such approximation, the 
FW and SHW fields are written as 

^{x,y,z)=A^ w (x,y)e' (k ^- mt \ 



(2) 



and 



E 2 (x,y,z) = B 0 (z)E 2 ,o(x,y)e' 



(fc2,oz — 2tof) 



(3) 



respectively. It is assumed that the FW field is non-depleted and has a 
constant amplitude A 0 , as in practice the fraction of FW power 
transformed to SHW power is negligible. Bo(z), which is the ampli- 
tude of the SHW field to be solved, is a slowly- varying quantity with 
respect to the interaction length z. E\#(x,y) and E,2fAx,y) are the 
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dimensionless transverse modal profiles of the lowest-attenuation 
FW and SHW waveguide modes with complex propagation wave 
vectors ki$ = k\ 0 + ik\ 0 and fc 2i o = k r 2 0 + ik' 2 0 , respectively. These 
two modal functions are normalized in practical operation with its 
maximum modulus set to unity. Substituting E 1; E 2 , and P <2) into 
Eq. (1), and using the well-established slowly-varying amplitude 
approximation, we arrive at the following relation 

dB 0 (z) 

2ik 2 fl^2,o{x,y) — - — 

dz ?- rl 



^- 5(x - x s ,y-y s )A 2 0 El 0 (x,y)e iAh 



(4) 



where Ak = 2k i0 — k 20 denotes the phase mismatch of the lowest- 
attenuation FW and SHW waveguide modes. The function 
8(x — x$,y— ys) means that only the inner wall surfaces of the air 
hole channel, whose positions are denoted by (xs,ys), make non-zero 
contributions to SHG signals. 

The method of moment is employed to solve Eq. (4). Both sides of 
Eq. (4) are multiplied by the transverse modal function ~E 2 Q (x,y) and 
then integrated in the transverse dimensions where the nonlinear 
susceptibility is non-zero. The following equation is obtained: 
dB 0 



dz 



-~iD(co)A z 0 Ze u 



(5) 



where D(cu) = (2uj 2 /} 2 ^)/(k 2 fit: 0 c 2 ), an d 3 is the physical quantity 
spatial overlap factor, which is defined as 
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Here S denotes the inner wall surface area of the air hole. The integ- 
rand is the product between the tensor polarization and the vector 
field of SHW. The normalized integration represents the degree of 
effective spatial overlap between the FW and SHW modal fields. 
Assuming the initial SHW field at z = 0 to be zero, we obtain the 
solution for Eq. (5) as 
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Having obtained the amplitude of SHW waveguide mode, the intens- 
ity of SHW at the exit surface (the back surface of the gold thin film, 
located at z = h, where h is the gold film thickness) is then given by 
Eq. (3) and Eq. (7). 
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The integration of Eq. (8) is upon the whole cross section area of the 
rectangular air hole. This average intensity is proportional to the 
observed far-field SHG signal radiating from the periodic array of 
subwavelength air hole array. 

Discussion 

According to Eq. (8), the SHG signal is determined by several factors. 
The first factor is the local field enhancement effect described by the 
amplitude A 0 of the FW mode Ei 0 (x,y). Note that this value does not 
directly correspond to the transmission power of FW, it just reflects 
the field within the hole channel. The second factor is the spatial 
overlap factor E, which is closely relevant with the modal profile 
similarity between the FW and SHW waveguide modes that take 
part in the nonlinear optical process of SHG. The third factor is 
the phase mismatch between the FW and SHW modes, which can 
also make a strong influence due to the strong dispersion nature of 
the metallic-waveguide mode, even if the thickness of the gold films is 
several times smaller than the wavelength. The fourth one is the 



attenuation strength of the SHW mode described by the quantity 
k' 2 0 . In essence, all these four factors can be closely dependent on the 
geometric shape of the air hole. The nonlinear coupled-wave theory 
has confirmed the validity of our above general description of the 
physical picture of SHG from these nonlinear nanophotonic struc- 
tures, and has quantified the contribution of each relevant physical 
mechanism. 

It is helpful to compare our theory with previous investigations in 
Ref. 31. In that work, the effective nonlinear response from the air- 
hole array is enhanced by more than one order of magnitude for a 
certain particular value of aspect ratio. This enhancement was merely 
attributed to the influence of the hole shape on the linear transmis- 
sion properties, in particular, the slow propagation of the fun- 
damental wavelength through the holes which occurs close to the 
hole cutoff. Other important characteristics on the nonlinear modal 
conversion between FW and SHW, as have been clearly analyzed in 
our theory, were not accounted for explicitly in Ref. 31. 

We first determined the characteristics of the lowest-attenuation 
FW and SHW waveguide mode, including their transverse modal 
profile function [Ei 0 (x,y) and E 20 (x,y)] and their dispersion, 
namely, the complex propagation wave vectors fc 1-0 and k 2 0 as a 
function of the air hole AR. This was done via the 3D FDTD numer- 
ical calculations. The dispersive nature of gold material was taken 
into account 36 . The simulation results are displayed in Fig. 3. In 
Fig. 3(a) and 3(b), the real and imaginary part of fc 10 
(ki$ = k r 10 + ik\ 0 ) and fc 2 ,o(fc 2i o = k 2 0 + ik' 2 0 ) are plotted for the 
800 nm FW and 400 nm SHW. For the FW, k\ 0 first keeps at a 
low level when AR increases from 0.5, then grows rapidly when it 
reaches around AR = 2.2. As k\ 0 describes the phase factor and is 
related to the phase index, the result might be a signature of the 
ignition of slow transport mode at AR = 2.2, which has been 
assumed in Ref. 31 as the major reason responsible for the shape 
resonance of SHG. However, k\ 0 keeps decreasing (almost linearly) 
when AR grows. Thus the variation feature of k\ „ does not match the 
picture of slow transport mode because in this dissipative material 
the ignition of slow mode usually means a greatly enhanced dissipa- 
tion of the mode. The simulation result of linear light dispersion 
properties, together with the measurement result of light pulse delay 
depicted in Fig. 2(b), strongly indicates that there is not obvious slow 
transport mode for FW in the rectangular hole with various AR. The 
dispersion curves of SHW, as depicted in Fig. 3(b) shows a very 
different behavior from that of FW. When the hole AR grows, k 2 0 
continuously increases without a distinct jump, while k' 2 0 decreases 
very rapidly. The results clearly show that the geometric shape, 
reflected in the aspect ratio, has a significant influence on the trans- 
port properties of both FW and SHW. The attenuation is much 
reduced for large AR than for small AR. This is not difficult to 
understand. As the incident FW is parallel to the y-axis, a larger 
AR means a larger x-axis space of the rectangular hole witnessed 
by the light, and therefore a smaller attenuation due to waveguide 
cutoff. 

Our numerical calculations show that the prominent field compo- 
nents of the waveguide mode Ei(x,y,z) andE 2 (x,jK,z)[seeEqs. (2) and 
(3)] under the excitation of the y-axis polarization incident light is 
the E y field, while other two field components are much smaller, 
similar to the incident light. As a result, for simplicity we may con- 
sider only the E y field to describe E\(x,y,z) and TL 2 (x,y,z). 
Figures 3(c)-(f) display the E^-field amplitude pattern of waveguide 
mode with wavelengths 800 nm and 400 nm at the entrance of the 
air hole for two specific hole shapes with AR = 2.3 and AR = 1.6. 
Note that a transverse mode has small variation of phase within the 
air hole area. The incident light has unity amplitude. These field 
patterns should correspond to Ei(x,y,0) and E 2 (x,y,0)and they 
involve the information of both the transverse mode profile func- 
tions [E 10 (x,y) and E 2 fi(x,y)] and the FW mode amplitude A 0 . It can 
be seen that the modal profile for the two hole shapes are similar for 
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Figure 3 | Physical properties of the FW and SHW waveguide modes, (a) Calculated imaginary part and the real part of the propagation constant at the 
wavelength of (a) 800 nm and (b) 400 nm, respectively, (c)-(f) £^-field amplitude pattern of at the wavelength of 800 nm [for the FW mode E l: o(x,y)] 
and 400 nm [for the SHW mode E2,o(x-,y)] for air hole with AR = 2.3 and 1.6. The dotted rectangular denotes the air hole boundary. The incident 
electric field is y-polarized and the dominant component of electric field is E y . All the field intensities are in unit of the incident field intensity. 



both the FW and SHW modes. Besides, the modal profile for FW and 
SHW does not differ much at the same AR value either, although FW 
and SHW seem to have slightly higher similarity for AR = 2.3 than 
for AR = 1.6. All the four modes have the maximum field intensity 
located at the center of the upper and lower x-axis edges of the hole. 
However, there is a prominent difference: The FW field amplitude is 
much stronger for AR = 2.3 than for AR = 1.6, indicating a much 
higher field enhancement of the FW mode for the AR = 2.3 hole. 

The above knowledge on the FW and SHW waveguide mode 
properties allows us to evaluate quantitatively the SHG signal intens- 
ity from various hole array samples. The calculation results based on 
Eq. (7) and Eq. (8) are summarized in Fig. 4 for a series of nano- 
structured gold thin films with various air hole AR. Figure 4(e) illus- 
trates the amplitude B 0 (h) of SHW waveguide mode at the exit plane 
of the thin film, which originated from the excitation by FW and 
evolution of SHW within the hole. Figure 4(f) illustrates the corres- 
ponding intensity of SHG signal emitting from the gold thin film. In 
these two figures, the results have been normalized to the value at AR 
= 0.5 (which is set to l).The calculated amplitude of SHW mode and 
SHG signal intensity all show a prominent resonant peak at AR = 
2.3. And there is a three-orders-of-magnitude variation of the SHG 
signal intensity for different AR samples. The theoretical result 
agrees well with our experimental data as shown in Fig. 2, as well 
as with those in Ref. 31. 



The theory developed in the last section has pointed out several 
factors that contribute to the overall SHG signals. The calculation 
results on the contribution from the following physical factors: (a) 
The amplitude of FW mode |A 0 |; (b) the spatial overlap factor |E|; (c) 
the FW-SHW phase mismatch factor j sin (A/ch/2)/(fc 2 , 0 Afc/2)|; and 
(d) the SHW attenuation factor exp( — k' 20 h) for various hole AR 
values are displayed in Figs. 4(a)-4(d), respectively. Figure 4(a) 
clearly shows that the amplitude of FW mode |Ao| almost grows 
monotonously (accompanied with some fluctuations) with the 
increasing AR of air hole, and the peak takes place at AR = 2.3 with 
the magnitude of 6.7 (normalized with the amplitude at AR = 0.5). It 
is clear from Eq. (7) and Eq. (8) that the SHW mode amplitude is 
B 0 cc |A 0 | 2 , while the SHG signal intensity is I$hg °c |Ao| 4 - As a result 
for this fourth power law, the modest enhancement on the excited 
FW mode amplitude \A 0 \ will lead to three orders of magnitude 
enhancement in SHG signal intensity. 

The FW-SHW spatial overlap factor |H| follows some complicated 
variation but overall it fluctuates around a high level of 0.7 for all the 
hole shapes (AR values). It can be said that this spatial overlap factor 
does not play an important role in shaping the SHG signal intensity 
for different hole shapes in the current nonlinear nanophotonic sys- 
tem. The reason is because that the modal profile of FW and SHW 
modes does not change much for the current air hole geometry. 
However, this factor might play a dominant role in other systems 
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Figure 4 | Theoretical results of SHG signals as a function of the air-hole aspect ratio. Panels (a)-(d) show the contribution from various physical 
factors to SHG: (a) The amplitude of FW mode A 0 ; (b) the spatial overlap factor HI; (c) the FW-SHW phase mismatch factor; and (d) the SHW 
attenuation factor. Panels (e) and (f) display the calculated amplitude of SHW mode and SHG intensity (linear scale, with the corresponding semi-log 
scale plot depicted in the inset), respectively, which show a prominent resonant peak at AR = 2.3 and a three orders of magnitude variation of the SHG 
intensity for different AR samples with a total area of S = 260 nm X 130 nm. 



where the FW and SHW modes have a very different modal profile 
and each of them depends sensitively on the geometry of the non- 
linear nanostructures. The above two factors shown in Fig. 4(a) and 
4(b) can be attributed to the transverse modal profile of FW mode 
and SHW mode, with \ A 0 \ being related to the coupling efficiency of 
FW mode with the incident light, while |S| being related to the 
conversion efficiency from FW to SHW. It can be said that these 
two factors reflect the local interaction property of FW with metallic 
nanostructures. 

In contrast, the other two factors shown in Fig. 4(c) and 4(d) 
are closely related with the dispersion nature of FW and SHW, 
namely with the subtle way how these two modes transport along 
the air hole channel. The FW-SHW phase mismatch factor 
j sin(Afch/2)/(/c 2 , 0 AA:/2)|[Fig. 4(c), all normalized to that at AR = 
0.5], which is well known to play a key role in conventional mac- 
roscopic nonlinear optics, e.g., in nonlinear crystals, is also very 
important in current nonlinear nanophotonic system. The phase 
mismatch factor almost decreases monotonously with increasing 
value of AR, with the value at AR = 2.3 about 2 times smaller than 
the value at AR = 0.5. The SHW attenuation factor exp( — k' 20 h) 
show an opposite variation behavior to that of the phase mismatch 
factor. It grows monotonously with increasing AR value, and gradu- 
ally reaches a saturation plateau when AR reaches 2.0. The value at 



AR = 2.3 is about 2 times larger than the value at AR = 0.5. 
Obviously the SHW attenuation factor also plays an important role 
in SHG from the air hole channel. It can be said that the above two 
factors reflect the global transport property of FW and SHW. 
Although they are all important for SHG, they counteract and cancel 
with each other, leaving the overall influence from the global trans- 
port property not a significant factor for the ultimate intensity of the 
observed SHG signal in the current nonlinear nanophotonic system. 
However, this might be an accidental event. In many other systems 
they altogether can be a significant factor in determining the intens- 
ity of SHG. 

There is some discrepancy between theory and experiment (see 
Fig. 2 and Fig. 4) in regard to the hole shape, i.e., AR value (2.3 in 
theory and 2.0 in experiment) at which the maximum SHG signal 
occurs. We attribute this discrepancy to the inconsistence in the 
geometric shape of practical samples with those adopted in numer- 
ical simulation. In simulation, we assumed perfect rectangular air 
holes with strictly vertical walls etched into a strictly perfect gold thin 
film. In practice, there exists considerable surface roughness (gold 
islands) at the inner and outer surfaces, the size of air hole cannot be 
perfectly matched with preset value due to the inaccuracy of FIB 
milling, there exists over-etching of gold in the front and back exit 
of the air hole, and the wall of air hole can be slightly inclined with 



SCIENTIFIC REPORTS | 3 : 2358 | DOI: 1 0. 1 038/srep02358 



6 



respect to the surface of gold thin film. All these deviations in the 
geometric shape of practical hole from the modeled one would mod- 
ify the waveguide mode properties and change the way how the SHG 
signal depends on the AR. The above assumptions have been con- 
firmed by our further numerical calculations for different values of 
rectangular hole area, which are illustrated in Fig. 5. Three values of 
air hole area are further considered, S = 260 nm X 140 nm, 270 nm 
X 140 nm, and 260 nm x 120 nm. Comparison with the situation of 
S = 260 nm X 130 nm (Fig. 4) shows that the maximum SHG 
intensity takes place at smaller AR for larger air hole area, and at 
larger AR for smaller air hole area. For S = 270 nm X 140 nm, the 
SHG resonance happens at AR = 2.0, which is very close to our 
experimental data shown in Fig. 2(a). In addition, the maximum 
SHG signal intensity also depends on the hole area. These results 
again strongly support the ideology that geometric shape plays a key 
role in SHG from nonlinear nanophotonic structures. 

Our above analysis has clarified the contribution from several 
relevant physical factors, in which the local FW field enhancement 
overall dominates the SHG signals, It is worthwhile to compare our 
findings with previous literatures, in particular Ref. 31. To account 
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Figure 5 | Calculated SHG intensity as a function of air-hole aspect ratio 
for different hole areas. The SHG intensity has been normalized to the 
value at AR = 0.6. 



for the beautiful experimental data of enhanced SHG, the slow light 
effect was employed, which leads to an enhanced effective nonlinear 
susceptibility of the hole array gold thin film that is adopted to 
describe phenomenologically the overall SHG response of the sam- 
ple. Although the slow light effect can contribute to local field 
enhancement, local field enhancement can be induced by many other 
mechanisms. As both the FW mode band diagram calculation and 
the pulse delay experiment have excluded the existence of a signifi- 
cant slow light effect, we assume that the higher coupling efficiency of 
external FW light energy into the FW waveguide mode plays the 
dominant role. In addition, the spatial overlap factor |S| and the 
SHW attenuation factor were not involved in the model considered 
in Ref. 31, and the FW-SHW phase mismatch factor was considered, 
but deemed to negligible due to the small film thickness compared 
with the wavelength. According to our analysis, these three factors 
each can make considerable effect on the SHG. 

In summary, we have theoretically and experimentally investi- 
gated the physical origin of SHG shape resonance effect in nonlinear 
nanophotonic structures made from periodic arrays of rectangular 
air hole with various aspect ratios perforated in gold thin film. We 
have developed a nonlinear coupled mode theory to account for the 
energy conversion from FW mode to SHG mode involved in the SHG 
process within the nanoscale channel. The theory is complete and 
accurate as it is based on Maxwell's equations for nonlinear nano- 
photonic structures. The complete theory reveals that several phys- 
ical mechanisms, including the FW mode excitation amplitude, the 
FW-SHW modal spatial overlap, and FW-SHW mode phase mis- 
match, and the SHW mode attenuation, can make significant con- 
tribution to the SHG signal intensity. They are all geometric shape 
sensitive and altogether act to induce the shape resonance effect of 
SHG from these nanostructured gold thin films. In short, the SHG 
shape resonance can be ascribed to the strongly shape dependent 
intrinsic phase and attenuation feature of the FW and SHW wave- 
guide mode and the external excitation efficiency of the FW mode. 

The theory agrees well with our experimental observations of 
shape dependent SHG and also with those reported in previous 
literatures that the SHG signal intensity can be three orders of mag- 
nitude larger at a specific hole shape than at other shapes. Both of our 
ultrafast pulse delay experiment and guided mode band diagram 
calculation exclude the long-held model that slow light effect would 
play a key role in such shape resonance of SHG. Our theory provides 
an accurate explanation for the long-emphasized but elusive shape 
effect in this puzzling problem. It may stimulate new angles and 
deeper insights to visualize general nonlinear nanophotonic pro- 
cesses and pave the way to engineering high-efficiency nonlinear 
nanophotonic structures. Our finding would inspire a very useful 
insight into other general nanoscale nonlinear optical problems that 
have been inconspicuous in conventional macroscopic nonlinear 
optics. 

Methods 

Sample preparation. Gold thin film of 160 nm thickness was deposited onto BK7 
glass substrates, with a 5 nm thick Ti adhesion layer, by means of magnetron 
sputtering method. The gold thin film is thus not a single crystal, but composed of 
poly-crystalline gold where domain boundaries exist on and right below all surfaces. 
A 20 X 20 square-lattice array of rectangular holes were milled on these films by 
means of focused ion beam milling method, with a spatial periodicity of 410 nm and a 
total area of8 um X 8 um (Fig. 1). In all the samples each hole is designed to have the 
same area. Due to the practical limitation of nanofabrication accuracy, the outer and 
inner surfaces of the nanostructured samples are not perfect atomic-level flat surface 
patterned with single-crystalline gold atom lattice, but involve many nanometer-scale 
roughened gold islands with different crystalline orientations (Fig. lc and 5e). Thus 
all surfaces of the film all possess the necessary spatial inversion symmetry breaking 
for SHG. In addition, all the surfaces have an amorphous morphology, leading to a 
uniform and isotropic second-order susceptibility for the surface gold atoms. 

Measurement of SHG signals. We used femtosecond laser pulse to investigate the 
linear and nonlinear light propagation. Laser (Coherent Inc.) pulses of 800 nm 
central wavelength, 80 MHz repetition rate, and 45 fs pulse duration were focused 
onto the sample with a f — 50.8 mm lens. The light spot on the sample had a diameter 
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of about 30 urn. The transmitted FW and SHW passed through a/ — 63.5 mm 
focusing lens. The samples were mounted onto a manual translation- rotation stage, 
with the y direction always parallel to the laser polarization. A flip mirror was 
employed to switch between the SHG and time delay measurements. For the SHG 
measurement a 20 mW average laser power was used, with the flip mirror in use. A 
low-pass edge filter with a transmission of T = 6 X 10~ 9 for the FW beam and T = 
98% for the SHG signal was used. 

Method to achieve ultrahigh temporal resolution. By averaging multiple scans, the 
temporal position of the cross -correlation peaks became very stable. The temporal 
resolving capability, or the experimental precision A, is explicitly determined by 

^precision = maxf<5 resolution, \J accuracy) 2 repeatability) 2 / N^j > where 5 are VaOOUS 

systematic errors, and N is the number of scanning. Our linear stage has a resolution 
of 0.02 um, accuracy of +/— 0.3 urn, and repeatability of +/— 0.1 um. The 
unevenness of the substrate thickness was also considered and measured carefully, 
which does not bring in a noticeable system error in time delay. Furthermore the 
method is insensitive to the pulse-to-pulse jittering and can survive the pulse duration 
broadening. These are curial for obtaining a decisive conclusion for such a precise 
experiment. Overall the total experimental error is about ±2 fs. 

Numerical simulations. The 3D-FDTD method was used for the calculation of 
electric fields. The period of the structure was kept as 410 nm. The source was set a.sy- 
polarized in all calculations. The wavelength of FW is 800 nm. We used periodic 
boundary conditions along x- and y-directions, and perfectly matched layers (PML) 
were added along z-direction. The grid size is set as 5 nm. The lowest- attenuation 
waveguide mode of FW and SHW was determined by allowing light transporting 
within a deep rectangular hole. The complex propagation wave vector (ki t Q and fc^o) 
could be retrieved by comparing the phase and amplitude between two appropriately 
selected adjacent cross-sectional planes {both sufficiently far away from the entrance 
of the hole) with a small separation distance, while the transverse modal profile 
function [E^ 0 (x,y) and E 2 ^(x,y)] could be obtained by recording the field profile in 
one plane. 
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